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Doped two orbital chains with strong Hund's rule couplings 
ferromagnetism, spin gap, singlet and triplet pairings 
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Different models for doping of two-orbital chains with mobile 5 = 1/2 fermions and strong, 
ferromagnetic (FM) Hund's rule couplings stabilizing the S = 1 spins are investigated by density 
matrix renormalization group (DMRG) methods. The competition between antiferromagnetic 
(AF) and FM order leads to a rich phase diagram with a narrow FM region for weak AF 
couplings and strongly enhanced triplet pairing correlations. Without a level difference between 
the orbitals, the spin gap persists upon doping, whereas gapless spin excitations are generated 
by interactions among itinerant polarons in the presence of a level difference. In the charge 
sector we find dominant singlet pairing correlations without a level difference, whereas upon the 
inclusion of a Coulomb repulsion between the orbitals or with a level difference, charge density 
wave (CDW) correlations decay slowest. The string correlation functions remain finite upon 
doping for all models. 

KEYWORDS: spin-1 systems, Mott-Hubbard transition, superconductivity, triplet superconductivity, ferromag 
netism, string order 
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A key problem in the research on the Mott-Hubbard 
transition consisits of the doping of a spin liquid ground- 
state with mobile holes. The antiferromagnetic (AF) 
5=1 Heisenberg (HB) spin chain is an example of 
such a spin liquid with a finite spin gap 1 ) and a hid- 
den topological order, the string correlation function. 2 ) 
Doping of 5 = 1 HB chains poses a number of interest- 
ing theoretical questions caused by the competition be- 
tween ferromagnetic (FM) order induced by the double 
exchange mechanism 3 ) and AF order which can result in 
completely different magnetic properties. Another rele- 
vant question is whether the spin gap is destroyed im- 
mediately upon doping or whether it persists, and the 
competition also determines which correlation function 
characterized by the single correlation exponent K p dom- 
inates in the thermodynamic limit of infinite system size. 

Only recently has it become possible to experimen- 
tally study mobile holes doped in a 5 — 1 chain in the 
system Y2- x Ca x BaNi05. 4 ) The two active Ni + orbitals 



are the 3d? 



and the almost localized 3cL 



bital. Strong Hund's rule couplings Jh between these 
two orbitals and very weak inter-chain couplings make 
this system an almost ideal 5=1 HB chain with a spin 
gap of A s ps 100A" in the undoped case. 5 - 1 By replacing 
off-chain Y 3+ with Ca 2+ , mobile holes can be introduced 
into the system and /iSR data show that these carriers 
indeed have spin 5 = 1/2. 6 ) The most interesting ex- 
perimental features upon doping consist of a reduction 
of the resistivity pdc by several orders of magnitude and 
that the temperature dependence of pdc n o longer can be 
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described by thermal activation over a barrier. 4 ) Further 
new states with spin 5 between 1 and 3/2 per impurity 
appear below the Haldane gap. 4 ' 

First theoretical studies on doped 5=1 chains con- 
centrated on localized 5=1/2 impurities. 7 ) For suffi- 
ciently weak couplings, they give rise to bound 5=1/2 
states below the Haldane-gap, 8 ) similar to the chain-end 
excitations. 9 ) However, it soon turned out that mobile 
carriers are needed for an accurate description of this sys- 
tem and the case of weak hopping 10 ) has been studied 
first. A more realistic effective Hamiltonian for fully mo- 
bile holes has been derived in 11 ) for one- and two-band 
models with infinitely strong Hund's rule couplings. A 
phase diagram which shows a FM region and phase sep- 
aration for large values of the coupling J between the 
5=1 spins has been obtained by exact diagonalization 
and (DMRG) studies of small systems. 12 ) This study 
also suggests a possible hole-bound region between the 
FM and phase separated ones. Based on DMRG calcula- 
tions for very large systems, we briefly have reported in 
a previous Letter on a new hierarchy of energy scales for 
a model with finite Hund's rule couplings including fully 
mobile holes and localized electrons in the lower band. 13 ) 
Ferromagnetic polarons formed by the double exchange 
mechanism are found to weakly interact in the lower en- 
ergy scale in the background of the gapped spin liquid 
formed at the higher energy. Moreover, doping-induced 
magnetization steps have been discovered for an exactly 
solvable model which includes biquadratic exchange. 14 ) 
The biquadratic terms are needed for the solution by the 
algebraic Bethe Ansatz. 

While all of the above mentioned studies are restricted 
to models with mobile electrons in the upper band only, 
the very different behavior of systems with and with- 
out a level difference between the two orbitals has first 
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been investigated by weak-coupling theories. 15, These 
studies show that the spin gap remains finite if there is 
no level difference between the orbitals, whereas the spin 
gap is destroyed immediately if the particle density be- 
tween the two orbitals differs. We have recently carried 
out the first study of this problem in the strong coupling 
regime by numerical methods and discovered dominant 
superconducting pairing correlations for systems without 
a level difference in contrast to dominant charge-density 
wave (CDW) correlations for systems with a level differ- 
ence. 17 -' This surprising result which exhibits completely 
different phases depending on the level difference empha- 
sizes the importance of a detailed study of realistic mod- 
els in the strong coupling regime. Unfortunately this is 
also the most difficult regime to investigate, and it is not 
clear whether additional terms need to be included in a 
realistic model. 

In this paper we present a detailed analysis of a realis- 
tic model in the strong coupling regime with strong but 
finite Hund's rule couplings Jp and fully mobile holes. 
Based on large scale numerical calculations by both the 
thermal DMRG (TDMRG) 18 ' and ground state DMRG 
method, 19,20 ) we will to cover what has been left out 
due to lack of space in our two previous publications. In 
addition, we include a model with a Coulomb repulsion 
and no level difference between the electrons in the two 
orbitals on the same site. 

We start with the discussion of a generic model for 
Y 2 - x Ca x BaNi05 which includes a level difference be- 
tween the two orbitals, and we will call this the asym- 
metric model in the following. In this case, the electrons 
in the lower 3d x 2_ y 2 orbital are almost localized, and all 
holes are doped into the higher lying 3d 3z 2_ r 2 orbital. 
We describe these mobile holes in the upper band by the 
following Hamiltonian: 

#kln = 'X 7, t'w-< ' - • "■<■■) v, (1.1) 

where the projection operator V prohibits doubly occu- 
pied sites, Cj i a is the particle creation operator on site j 
in the orbital i with spin a, and the index i = 1 denotes 
the upper orbital and i = 2 the lower orbital. Due to vir- 
tual hopping processes there exist AF couplings J > 
between nearest-neighbors in the upper band, and we in- 
clude further AF couplings Jd > between the electrons 
in different orbitals on neighboring sites as follows: 

H ri = J Y1 (Sj-iSj+ht ~ \". l,»J (1-2) 

#diag = Jd"^ ^,1^7+1,2 + %2<Sj+l,l ~ 2 n J' in J+ 1 > 2 ^ ' 

(1.3) 

where the indices are the same as before, njj = c\ i i-Cj »f+ 

j i i c J>«»i' ano - * ne rest °f * ne no t a ti on is standard. The 
strong Hund's rule coupling Jh < between the elec- 
trons in different orbitals on the same site reads: 

^fm = Jh ^ Sj,iSj t 2, (1.4) 
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and is responsible for the formation of the S = 1 spins. 
Finally the Hamiltonian for a model with a level differ- 
ence is expressed as: 

# as = H FM + H dmg + H$ + H^l (1.5) 

This model is easy to understand in two limiting cases: 
first at half filling, where it can be mapped to a Haldane 
S = 1 chain with a spin liquid ground state and a finite 
spin gap. The second case is away from half-filling at 
J = Jd = 0, where all the spins are ferromagnetically 
aligned due to the double-exchange mechanism in order 
to gain kinetic energy. Hence, for finite values of J, Jd 
there is a competition between FM order induced by the 
double exchange mechanism and the spin liquid ground 
state. For \Jh\ J,Jd each hole is surrounded by a 
small FM cloud created by the double exchange mecha- 
nism, and we will call this a polaron in the following. An 
intriguing question on the doped system for small J and 
Jd is the possiblity of triplet pairing. In this paper we 
show a strong enhancement of the triplet pairing corrlcla- 
tions. However, for J, Jd > O.lt, ferromagnetic polaron 
becomes only a weak perturbation of the underlying spin 
liquid which remains intact. Among these polarons 2k p 
and 4fcp charge density wave (CDW) order is then sta- 
bilized. A hierarchy of energy scales is established in 
the spin sector with the smaller scale given by the gap- 
less, lowest-lying interactions among the polarons; and 
the second, larger energy-scale consists of the spin liquid 
background. 

On the other hand, we can think of a model with a 
very small or no level difference at all. This situation is 
more likely in a ladder-system, and we will refer to this 
model as the symmetric model in the following. From the 
symmetric model we include hopping and AF couplings 
J > in both orbitals, and the Hamiltonian is given by: 

H sym = H FM + H diag + £ + H®) . (1.6) 

i 

In this model the double exchange mechanism is ex- 
pected to be less effective, since hole pairs on the same 
rung can gain the large Hund's rule coupling Jh- How- 
ever, in many realistic situations there might be an ad- 
ditional Coulomb repulsion between the electrons on the 
same rung, and we will include this by the following ad- 
ditional term: 

H rcp = U y^l> (1-7) 

3 

This Coulomb repulsion will reestablish the competition 
between the double-exchange mechanism and the spin 
liquid ground state, while preserving the same particle 
density in both orbitals. 

For all the models we restrict our investigations to the 
case of strong Hund's rule couplings \Jh\ S> |J|, \ Jd\, \t\, 
and if not otherwise mentioned we always set —Jh = 
lOt = 20 J = 20 J d , and U = \J H \ for the model with 
a Coulomb repulsion between the electrons on the same 
site. 

The rest of this paper is organized as follows: we start 
with a discussion of the numerical methods first. In the 
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next section we investigate various physical quantities in 
detail, and we begin with the discussion of the spin gap. 
There we see that the spin gap is destroyed immediately 
upon doping for the asymmetric model, while it remains 
finite for the symmetric models. Consistent with these 
findings we obtain finite values of the magnetic suscep- 
tibility x f° r T — > for the asymmetric model. We 
continue with the study of several correlation functions 
obtained by the DMRG method, and from the correla- 
tion exponent K p obtained by fitting to the asymptotic 
form of a Tomonaga-Luttingcr liquid in the gapless and a 
Luther-Emery liquid in the gapful case, we find that su- 
perconducting pairing correlations dominate in the ther- 
modynamic limit of infinite system size for the symmetric 
model without Coulomb repulsion, while CDW correla- 
tions dominate in the two other cases. Finally, we in- 
vestigate ferromagnetism and triplet pairing correlations 
for small values of the AF couplings J, J d , and the last 
section of this report is devoted to the conclusion. 

§2. Numerical Methods 

In order to take the strong correlation effects into ac- 
count properly, unbiased and approximation free meth- 
ods are required, and large system-sizes are needed to 
avoid finite-size effects. For the investigation of the ther- 
modynamic properties, we use the novel TDMRG me- 
thod. 18 - ) This method is based on the combination of 
the transfer matrix method 21 ) and the DMRG method. 
The TDMRG method allows us to reach temperatures as 
low as T = 0.02i in the thermodynamic limit of infinite 
system size; and the only two sources of errors stem from 
the truncation errors of the DMRG procedure and the 
finite size of the Trotter time steps. 

In order to apply the transfer matrix method we first 
map a one-dimensional quantum system to a two-di- 
mensional classical system by the Trotter-Suzuki decom- 
position. 22 ) For this purpose we decompose the Hamil- 
tonian into two parts H = JJ dd + H cvcn where each 
-ffodd(cvcn) is a sum of commuting terms 7J dd = 
YhII h2i-i,2i and H cvcn = Ya=\ h-2i,2i+i- We then write 
the partition function of a system with L sites as 

Z L = tie- (iH = lim Z N . L (2.1) 

N—>oo 

with Z N x = tr(e-' 3 ^d/A r e -/3Hove„/JV^ Ar and ^ = x j T 
With the definition of the virtual transfer-matrix % — 

( e -0fc 2 *-l,2*/JV e -/3/»2 ( ,2*+l/iV) JV wc l = 

lim n oo^ Y\ i= i L %- By applying the Trotter-Suzuki 
interchangeability theorem 22 ) we can exchange the limit 
of the system size with the limit of the Trotter number 
and obtain the free energy density / in the thermody- 
namic limit of infinite system size 

/ = UmN^oolimL^oo (—^fj mtr ^ 

^ * / z=l,L 

= UrriN^oo (^--jjj In^max, (2.2) 

where A ma x is the maximum eigenvalue of %. We can 
then calculate the maximum eigenvalue A max numerically 
and obtain physical quantities from derivatives of /. In 



the actual calculations, several thermodynamic quanti- 
ties such as the internal energy u or the magnetization 
m are obtained directly from the maximal eigenvector 
of 7i. 21,23 ) For a fixed value of the Trotter time steps 
At = Pq/Nq, we then lower the temperature T by in- 
serting additional time steps along the imaginary time 
direction between the system and environment parts of 
the transfer matrix %. At the same time the dimension 
of the transfer matrix % is kept tractable by applying 
the DMRG method and truncating the states with the 
smallest weight in the density matrix. The dimension of 
the Hilbert space can further be reduced by using sym- 
metries such as the spin conservation symmetry and per- 
forming the calculations in the subspace of zero winding 
number. 

In our calculations we usually keep M = 80 states for 
each of the system and environment blocks of the DMRG 
procedure, and estimate the truncation errors by com- 
paring to the results with M = 40 and 60 states. For 
M = 60 states, the truncation error e is always smaller 
than 10~ 3 . Another source of errors stems from the fi- 
nite size of the Trotter time steps At. Basically one 
can extrapolate to At — > by performing calculations 
for several different values of At. However, for doped 
systems, the TDMRG calculations need to be performed 
in the grand canonical ensemble, and the chemical po- 
tential \i is fixed for each DMRG iteration. Thus wc 
need to extrapolate to constant particle density n af- 
terwards, and for a large number of Trotter time steps 
N, the errors introduced by this extrapolation to con- 
stant particle density are in general larger than the er- 
rors caused by the finite size of the Trotter time steps. 
If not otherwise mentioned, we have therefore performed 
the calculations for fixed At = 0.2i and neglected the 
errors caused by finite Trotter time steps. We have also 
used a re-biorthogonalization algorithm if we encounter 
numerical instabilities in the TDMRG procedure. 24 ' 

Complementary to the TDMRG method, we have used 
the ground state DMRG method to investigate the cor- 
relation functions and ground-state energies. The main 
advantage of this method is that it allows very accurate 
calculations of these quantities for large ID systems. We 
have investigated systems of up to 2 x 256 sites and kept 
up to M = 1400 states for the system and environment 
block each. The truncation errors encountered in our 
calculations are usually smaller than e < 10~ 8 . All cal- 
culations have been performed on Alpha workstations at 
the ISSP. 

§3. Results 

3. 1 Spin gap and magnetic properties 

We start our investigations with the spin gap proper- 
ties, and determine whether it remains finite upon dop- 
ing for the different models. In the undoped case, both 
models can be mapped to the Haldane S = 1 chain for 
Jh 3> J, Jd- The effective coupling J e g between the S = 
1 spins is J^fi = 3 J/4 for the asymmetric and J^g m = J 
for the symmetric model by second order perturbation 
theory, thus we expect a spin gap of A s = 0.41050(1) J e ff 
in both cases. 

By finite-size scaling we have numerically determined 
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£sym 


^sym+rc 


^sym+rc 
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Xo 





0.205(5)* 


6.05(5) 


0.205(5)* 


6.05(5) 


0.1504(5)* 


6.05(5) 







0.0625 


0.148(5)* 


6.1(2) 


0.03(1)* 


15.4(5) 





7.9(2) 


0.244 


2.61 


0.125 


0.120(5)* 


6.7(2) 


0.02(1)* 


20.5(5) 





11.2(2) 


0.201 


3.17 



Table I. Spin gap obtained by finite-size scaling analysis for the symmetric model (left), symmetric model with Coulomb repulsion 
U = 10* (middle), and the asymmetric model (right) for — Jh = 10* = 20J = 20Jd- Also listed is the spin-spin correlation length £, 
and for the gaplcss asymmetric model the spin velocity v a and magnetic susceptibility xo at T = 0. The definition of £ is not simple 
for the cases with A s = 0, see text for details. 
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Fig. 1. Finite-size scaling of the spin gap A s for the different 
models with -J H = 10* = 20 J = 20 J d , and U = 10* for the 
model with Coulomb repulsion, and different hole densities of 
= 0, rife = 1/8, and «/, = 1/16. 



the spin gap from A s = liniT^o A S (L; N — Ln), where 
A S (L;N = Ln) = A S {L;N) = E {L;N;S Z = 1) - 
Eo(L; N; S z = 0) with E (L; N; S z ) denoting the ground- 
state energy of the system with N particles on L sites and 
total spin component along the z-direction S z . The re- 
sult of the finite-size scaling analysis by DMRG calcula- 
tions (see Fig. 1) is in excellent agreement with the above 
mapping, and we obtain A s = 0.41(1) w 0.154(5)t in 
the asymmetric case and A s = 0.41(1) Jgg™ ~ 0.205(5)t 
for the symmetric model. In agreement with refs. 15 ' 16 -* 
we also find completely different behavior for both mod- 
els upon doping, as is evident from Fig. 1 and we have 
listed the results in Table I. For the asymmetric model, 
the spin gap is completely destroyed already for hole den- 
sities of nh — 0.0625 and rih = 0.125, where nh denotes 
the hole density in the conduction band. On the other 
hand, the gap is reduced, but remains finite for the sym- 
metric model for the same hole densities nh — 0.0625 and 
nh = 0.125, where we set nh = 1 — n with n denoting the 
number of electrons per orbital for the symmetric mod- 
els. The reduction of the spin gap is stronger in the sym- 
metric case including a Coulomb repulsion U = Jh, but 
appears to remain finite. However, great care needs to be 
taken with the boundary conditions of the DMRG cal- 



culations in order to avoid the degenerate 5 = 1/2 spin 
excitations at the end of the open chain, which would 
result in spin gaps exponentially small with the system- 
size L. 19 ) By enumerating the sites of the open chain 
from i = 1 to L, we have used the boundary condition 
which enforces a S — 1/2 spin at site i = 1 and S = 1 at 
site i = 2, and similarly at the other end of the chain. If 
we do not enforce a S = 1 spin next to the S =1/2 spins 
we find that holes are trapped at the boundary and the 
finite-size scaling becomes unreliable or impossible for 
smaller system sizes because of chain-end excitations. 

From the finite-size scaling we can also determine the 
spin velocity and magnetic susceptibility \ at T = of 
the gapless asymmetric model by noting that the small- 
est momentum in the finite open chain is fc m j n = ir/L. 
The spin velocity v a can then be obtained from A S (L) = 
Vakmin, and the magnetic susceptibility is given by xo = 
where K a = 1 due to the SU(2) symmetry. We see 
from Table I that the magnetic susceptibility \ 1S rather 
large already for small values of hole doping. 

Next we compare these findings for the magnetic sus- 
ceptibility with the TDMRG results shown in Fig. 2 for 
the asymmetric model for — Jh = lOt = 20J = 20 Jd, and 
we refer to Fig. 1 of ref. 17 for the symmetric model. For 
the symmetric model we have set — Jh — > oo to be able 
to keep enough states in the TDMRG algorithm (with 
M = 100 states). In the inset of Fig. 2 we show the 
convergence of the TDMRG method for several choices 
of the Trotter time steps At and numbers of states kept 
for the example of nh = 0.125. Keeping only M = 45 
states is clearly not sufficient, and this curve deviates 
from the other curves already at rather high tempera- 
tures around T w O.St. The other curves with M = 60 
and different values of At show good convergence, but 
at the lowest temperatures an accurate estimate of the 
errors gets more difficult. In the figure for all hole den- 
sities we therefore only show \ down to temperatures 
where we estimate the truncation errors to be smaller 
than 1% and interpolation errors to constant particle 
density also smaller than 1%. The first observation 
is that both the symmetric and asymmetric model show 
a strong enhancement of \ at temperatures near the gap 
value of the undoped chain T oc J/2. This strong en- 
hancement of x lS caused by the creation of a small 
FM cloud around each mobile hole by the double ex- 
change mechanism, as the holes can gain kinetic energy 
by FM alignment of the neighboring spins. At even lower 
temperatures, the magnetic susceptibility is suppressed 
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Fig. 2. Magnetic susceptibility x of the asymmetric model for 
-J H = Wt = 20J = 20J d . The T = values arc obtained from 
the spin velocities obtained by a finite-size scaling analysis. In 
the inset we show the convergency of the TDMRG method for 
different sizes of the Trotter time steps At and different numbers 
of states kept M for «/, = 0.125. 



for the symmetric model due to the spin gap and we 
find x - * for T — > 0. In contrast, \ IS finite and 
rather large at T = for the asymmetric model, indi- 
cating the formation of larger FM moments due to the 
double exchange mechanism (see ref. 13). However, the 
ground state is spin singlet for both models, as we have 
tested with the DMRG method by the calculation of 
a(i,j) = SfSI - l/2(S+Sj) which vanishes for a ro- 
tationally invariant ground state and the result is zero 
within the numerical precision of the DMRG for this 
quantity for both models (\a(i,j)\ < 10~ 6 Vz, j). 

3. 2 Spin correlation junction 

Complementary information about the magnetic prop- 
erties can be gained from the spin-spin correlations S sp (i— 
j) = (SfSj). For gapful systems we expect the spin-spin 
correlations to decay exponentially, whereas we expect 
power-law decay for gapless systems. By the DMRG 
method we have measured the spin-spin correlations for 
the different models on very large chains of up to 2 x 
256 sites, and we show some of those calculations in 
Fig. 3. For the undoped system, the spin-spin correla- 
tions clearly exhibit the expected exponential decay, and 
by fitting to the asymptotic form S sp (l) oc (— l)'Z _1 / 2 e~'/£ 
we obtain £ = 6.05(5) in agreement with the results for 
the Haldane S = 1 chain. 19 ) For both gapful sym- 
metric models we expect exponentially decaying spin- 
spin correlations also in the doped case. However, in 
Fig. 3 we observe an additional spiral-order in S?Sj on 
top of the exponential decay upon doping, and we fit 
to S ap (l) oc cos(2kpl)e^ 1 ^ in that case. The correlation 
lengths listed in Table I demonstrate that correlation 
lengths £ w 6.1(2) for n h = 0.0625 and £ « 6.7(2) for 
nh = 0.0625 are very close to the undoped case in the 
absence of the Coulomb repulsion, while the correlation 
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Fig. 3. Spin-spin correlation functions 5? 5? at low hole doping 
nh = 0.125 for the different models calculated by the DMRG 
method on 2 X 128 site systems and —J H = 10t = 20 J = 20 J<j, 
and U = —Jh f° r the symmetric model with Coulomb repulsion. 
For comparison we also show the undoped symmetric system. 



length £ increases by a factor of 2.56 to £ w 15.4(5) for 
n h = 0.0625 and 3.41 to £ « 20.5(5) for n h = 0.125 
at U = lOt. Rather surprisingly, also for the gapless 
asymmetric model the spin-spin correlations seem to de- 
cay exponentially with correlation lengths £ = 7.9(2) for 
n h = 0.0625 and £ = 11.2(2) for n h = 0.125. This ex- 
ponential decay of the spin-spin correlations at interme- 
diate distances reflects the underlying spin-liquid back- 
ground, and we expect a crossover to power-law decay 
at larger distances. However, since this length scale is 
given by the polaron-polaron distance, we can not dis- 
criminate between power-law and exponential decay due 
to the fact that even our largest systems contain few po- 
larons only, and much larger systems would be required 
for that purpose. 

This hierarchy of energy scales in the spin-sector of the 
asymmetric model should also be visible by two peaks in 
the specific heat cy. First, we expect a large, broad 
peak around the gap-value A s of the undoped S = 1 
chain stemming from the spin-liquid background, and 
a second peak at low temperatures reflecting the lower 
gapless AF interactions among the polarons. This two- 
peak structure is evident in Fig. 4 of ref. 13 by a broad 
peak around T w J reflecting the spin-liquid background 
and from the steep increase of cy signaled at the lowest 
temperatures. We estimate an energy scale of T < 0.02t 
for the magnetic interactions among the polarons. 

3. 3 String correlation function 

The hidden Z 2 x Z 2 symmetry of the Haldane S = 
1 chain is revealed by the string correlation function 
9(1) = ((S*j (nllt^\e™ s Z)(S* o+l )), which quickly 
approaches the value g(l) w —0.374 for I ^> 1 in the 
undoped case. 2 ) As a natural extension for the doped 
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Fig. 4. String correlation function g(xQ,x) = ((X]j=i 2 &x i) 

(n^ 0+ i,,=i,2^^) £1=1,2^,1)) ^ low hole doping for 
the different models with — Jh = 10t = 20J = 20Jd and (7 = 
— Jh for the symmetric model with Coulomb repulsion obtained 
by DMRG calculations on systems with 2 X 128 sites. 
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Fig. 5. Friedel oscillations in the charge density distribution due 
to the open boundary conditions for the different models with 
-J H = 10t = 20J = 20J d at low hole doping n h = 0.125. In the 
inset we show the corresponding Fourier transform. 



models we consider 

g{x) = / ( £ S' XOti ) I J] ) ( ]T ■ 

\ *=1,2 \fc=x +l,j=l,2 / 1 = 1,2 I 

(3.1) 

and the results are shown in Fig. 4. In the undoped case 
we quickly recover the original result of g(l) —0.374 
for I » 1. The value of g(l) is only slightly reduced to 
g(l) w -0.31(2) for n h = 0.0625, and g(l) w -0.26(2) for 
nh = 0.125 for the symmetric model without Coulomb 
repulsion in agreement with the weak-coupling result of 
rcf. 15. Combined with the very small deviation in the 
spin correlation length £ from the undoped case, we con- 
lude that the spin structure of the symmetric model is 
very similar to that of the undoped S = 1 Haldane chain. 
The suppression of the string correlation function g(l) 
is much stronger if we include the Coulomb repulsion 
or for the asymmetric model and we obtain g(l) < 0.1 
for both of these models for nh = 0.125. In addition, 
the latter two models also exhibit large oscillations in 
the amplitude of g(x), caused by the Friedel oscillations 
in the charge density due to the open boundary condi- 
tions. By comparison g(l) is constant within 1.2% for the 
symmetric model without Coulomb repulsion apart from 
boundary effects. In fact, for the asymmetric model we 
again expect a crossover to power-law decay of g(l) after 
sufficiently many oscillations. 

3.4 Friedel oscillations and correlation exponent 

Friedel oscillations are induced in the charge density 
due to the open boundary conditions. For gapful Luther- 
Emery and gapless Tomonaga-Luttinger liquids, we can 
then obtain from these oscillations the correlation ex- 
ponent K p which determines the dominant correlation 
functions in the thermodynamic limit. 

We first consider the Fourier-transform of the charge 
density shown in the inset of Fig. 5. For both of the 



symmetric models we find only one peak at k = 2nir. 
Since the lower band is fully occupied for low hole den- 
sities, there are iV^ p = 2(1 — n)L holes in the upper 
band and its particle density is n up = 2n — 1. For 
large Fermi volume, the Fermi vector is thus given by 
2kp — 7r + (2n — 1)tt = 2nn and the above peak is 
compatible with 2kp oscillations. For the asymmet- 
ric model we find only one peak at k = n^ix at the 
smallest hole density n h = 0.0625, but two peaks at 
k = n h ir and k = 2n h n for n h = 0.125 as is evident 
in Fig. 5. Since we denote the hole density in the con- 
duction band with nh in this case, we obtain a Fermi 
vector of 2fc_p = n + (1 — nh)ir = nhir for a large Fermi 
volume and the above peaks are compatible with 2fc_F 
and 4&f oscillations. 

After having identified the 2kp and 4fci? fluctuations, 
we determine the correlation exponent by fitting to the 
Friedel oscillations of an impurity potential. In the gap- 
ful case they are given by S n (x) oc C\ cos(2k F x)x~ Kf ' + 
C 2 cos(4k F x)x- 2K " and S n (x) oc d cos{2k F x)x- {1+K "^ 2 + 
Ci cos(AkFx)x~ 2K f> in the gapless case. 25 ) For both gap- 
ful and gapless universality classes, a value of K p > 1 
leads to dominant pairing correlations in the thermody- 
namic limit, while the CDW correlations decay slowest 
for K p < 1. From Fig. 5 it is evident that the ampli- 
tudes of the symmetric model without Coulomb repul- 
sion decay much faster than for the other models where 
they are almost constant. For the fit we have to dis- 
card the states next to the boundary because of trapped 
states. For the symmetric model without Coulomb re- 
pulsion we find by fitting to the gapful Luther-Emery 
liquid K p w 1.5 ± 0.05 for = 0.125, showing dom- 
inant pairing correlations, and K p w 0.41 ± 0.08 upon 
the inclusion of a Coulomb repulsion of U = lOt, giving 
dominant CDW correlations. For the asymmetric model 
we fit to the gapless Tomonaga-Luttinger liquid and find 
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Fig. 6. Singlet pairing correlations P;(/)P/ +X (/') with P> (f) = 
7j^( c !,T c !+/,| -c !,i c !+/,T) for ( a ) symmetric model, (b) symmet- 
ric model with Coulomb repulsion U = lOt, and (c) asymmetric 
model at n h = 1/8 and -J H = 10* = 20J = 20J d . 



Fig. 7. Triplet pairing correlations Pi(f)P< + J.f) with P/(/) = 
-i=(c| ^c\ + f |+ c l x c I+/ ^ or symmetric model, (b) symmet- 
ric model with Coulomb repulsion U = lOt, and (c) asymmetric 
model at n h = 1/8 and -J H = lOt = 20 J = 20 J d . 



iip w 0.51 ± 0.05, also at nh = 0.125. Pair formation 
can be explained by the gain of the large Hund's rule 
coupling for holes pairs formed on the rungs of the sym- 
metric model. This energy gain is not possible for both of 
the other models, and in contrast the holes may gain ki- 
netic energy by the double exchange mechanism in these 
cases. In the low hole doping region we can determine 
hole pairing from the pair binding energy obtained from 
Ap a i r = 2Ei — E — E 2 , with E n denoting the ground- 
state energy with n holes. For the symmetric model 
without Coulomb repulsion we confirm hole pairing by 
a large binding energy of A pa ; r « 2.29i. Including the 
Coulomb repulsion, the holes repel each other and we ob- 
tain Ap a i r w — 0.07i, while the asymmetric model shows 
a very small positive binding energy of A pa ; r « 0.016t 
near half-filling. Pair binding has also been found in 
ref. 12 ' in a similar model. 

3.5 Pairing correlations 

Independently of the Friedel oscillations in the charge 
density, we can also determine the correlation exponent 
from the pairing correlations defined as Pi(f')P} +x (f), 

where P}(f) = -^{c\^. + ^ T C !,i c !+/,T^ 6ith6r denotes 
the singlet (— ) or the triplet pair (+) creation operator. 
The pairing correlations also give further information on 
the form-factor /. For gapful models belonging to the 
universality class of Luther-Emery liquids, the pairing 
correlations decay as Pi{f)P} +x {f) ai _1 ^ p , 26 ' whereas 
for gapless Tomonaga-Luttinger liquids the asymptotic 
form is given by Pi{f)P} +x {f) K ln(x)- 1 - 5 x" 1 - 1 /' R ''' + 
A 2 coa(2k F x)x- K "- 1 / K " 2 ^ 

In Fig. 6 we show the singlet pairing correlations ob- 
tained from our DMRG calculations on systems with 
2 x 128 sites. From the numerous possibilities of the 
form-factor / for the symmetric model, we restrict the 
discussion to those with the largest amplitudes in the fol- 



lowing. In excellent agreement with the previous results 
for K p we find K p s=s 1.55 ±0.05 for the symmetric model 
without Coulomb repulsion at nu = 0.125. The form fac- 
tor of the pairs with largest amplitudes / = (1, 1), (2, 1) 
and (3, 1) are consistent with a d x _ y symmetry analogue 
for a ladder. Surprisingly there are also almost no 2Uf 
oscillations in the pairing-correlations for the largest am- 
plitudes, despite the Friedel oscillations in the charge 
density. 

Upon the inclusion of the Coulomb repulsion, the pair- 
ing correlations decay much faster and in agreement with 
the previous result from the Friedel oscillations we ob- 
tain K p w 0.42. Finally by simultaneously fitting to 
/ = 2,4,6, and 8 for a Tomonaga-Luttinger liquid, we 
obtain K p « 0.51 ± 0.05 for the asymmetric model at 
nh = 0.125. The form factor of the pairing correlations 
with the largest amplitudes is / = 8 and / = 6 in this 
case. Rather extended pair formation is also evident 
from a corresponding, weak substructure at a distance 
of d = 8 in the hole pockets of the charge density distri- 
bution in Fig. 5. 

The triplet pairing correlations decay exponentially 
for both universality classes. The correlations for the 
form factors / with the largest amplitudes are shown 
in Fig. 7 for systems with 2 x 128 sites and P}(f) = 
-j= (A . A - + ct . A ) . The correlation functions for 

V2 *> I i+f,l J >4- i+/,T 

the other triplet states are similar. At very short dis- 
tances d < 5, triplet pairing correlations for rung pairs 
with / = (0,1) and / = (1,1) have the largest ampli- 
tudes for the symmetric model without Coulomb repul- 
sions, but since they decay exponentially fast, the ampli- 
tudes for singlet pair correlations become larger for dis- 
tances d > 5. Including the Coulomb repulsion U = lOt, 
the singlet correlations always have larger amplitudes. 
For the asymmetric model, the triplet correlations have 
the largest amplitudes for short distances d < 5 with 
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a form-factor of / = 4, but again, the singlet correla- 
tions become dominant at larger distances d > 5 due the 
exponential decay of the triplet correlations. 

3. 6 Ferromagnetism and triplet pairing 

Ferromagnetism induced by the double-exchange mech- 
anism is to be expected for the asymmetric model and the 
symmetric model with Coulomb repulsion at J = Jd = 0, 
and might extend to small, but finite values of J, Jd- 
Rather large FM regions have been identified in the phase 
diagrams of refs. 12,28 for related models. 

In a first approach we determine the FM region by 
comparing ground-state energies in different total S* ot 
subspaces calculated by the DMRG method. In the FM 
phase, all the spins are fully aligned and the ground state 
is compatible with S^ ot = (2L — Nh)/2, while for singlet 
ground state with Nh even, S* ot = is lowest in energy, 
where Nh denotes the number of holes in the system 
of length L. Because of the number of subspaces with 
different S^ ot becomes very large for large system sizes, 
we can do such calculations only for rather small sys- 
tems, and we have used L — 32 for our calculations. For 
the asymmetric and the symmetric model with U = 5t, 
S* ot = is lowest in energy for J = J d > O.lt, and S* ot = 
(2L — Nh)/2 is lowest for J = Jd = 0. For intermediate 
values < J, J d < O.lt some < S^ t < (2L - N h )/2 
is lowest in energy. While these results clearly demon- 
strate a singlet groundstate for J = Jd > O.lt and a FM 
groundstate for J = Jd = 0, the results are not con- 
clusive for the intermediate region, due to the special 
boundary conditions required to suppress the chain-end 
excitations. 

An alternative approach to determine the FM region 
is possible by the spin-spin correlations. Working in 
the S* ot — subspace, the spins are ferromagnetically 
aligned in the x — y plane in the FM phase. In Fig. 8 
we show the S + {2)S~(x) correlations for the asymmetric 
model and the symmetric model with U = 5t calculated 
by the DMRG method for a system size of L = 32 sites. 
In the asymmetric model, all the spins are fully ferro- 
magnetically aligned for J = Jd = 0. Ferromagnetic 
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alignment of the spins is still evident for J = Jd = O.Olt, 
however at the end of the chain, the last spin is anti- 
ferromagnetically aligned. For J = Jd = O.lt, we find 
AF alignment compatible with the singlet groundstate. 
The same situation is encountered also for the symmetric 
model with U = 5t, where we find fully polarized spins 
for J = Jd = O.Olt, a few spins at the boundary with 
AF aligment for J = Jd = 0.05i, and AF alignment for 
J = Jd = O.lt, where the period of the AF alignment 
corresponds to 2Uf oscillations. Without Coulomb re- 
pulsion, the FM phase seems to be slightly reduced, since 
the strong pair binding on the rungs reduces the double 
exchange mechanism. Already for J = Jd = O.Olt a few 
spins are antiferromagnetically aligned in this case, but 
the bulk is still ferromagnetically aligned. The AF align- 
ment of the spins at the boundaries is expected to be a 
finite-size effect only and can be explained by the de- 
creased mobility of the holes near the boundaries, which 
in turn reduces the double-exchange mechanism and fa- 
vors AF order. At Jh = 10t we therefore expect a narrow 
FM region for all models for < J = J d < O.lt. 

The possiblity of triplet superconductivity in the FM 
phase is a question of great interest and we have inves- 
tigated this problem by the pairing correlation shown in 
Fig. 9. As is evident from the figure, the triplet pairing 
correlations for the symmetric model with Jh = —10t, 
and J = Jd = O.Olt at nh = 0.125 for a triplet pair 

P i(?) = ^(4t c !+/U +c li c I+/ P formcd 011 a mn S arc 
roughly three orders of magnitude larger than the largest 

singlet pairing pairing corrletions, and at least for this 
system size, do not seem to decay faster than the sin- 
glet pairing correlations. Also for the asymmetric with 
the same parameters, the triplet pairing correlations are 
roughly three orders of magnitude larger than the singlet 
pairing correlations, and do not seem to decay faster than 
the singlet pairing correlations. From the system sizes in 
our study it is not possible to distinguish between expo- 
nential or power-law decay, but at least for the 2 x 32 
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systems, dominant triplet pairing seems possible. Larger 
system sizes are needed in order to determine which cor- 
relation function decays slowest in the thermodynamic 
limit of infinite system size. 

§4. Conclusions 

In our study we have investigated three realistic mod- 
els of AF 5=1 Heisenberg chains doped with mobile 
5 = 1/2 fermions in the strong coupling regime: a model 
with a level difference between the two orbitals form- 
ing the 5=1 spins, a model without level difference, 
and a model without level difference including a strong 
Coulomb repulsion between the electrons in the two or- 
bitals on the same site. Our investigation by numerical 
methods shows very different physical properties of these 
models. 

For the asymmetric model with a level difference be- 
tween the two orbitals, the spin gap is destroyed im- 
mediately upon doping, while it remains finite for both 
symmetric models with equal particle density in both or- 
bitals. This collapse of the spin gap in the asymmetric 
model is caused by AF interactions among the polarons. 
The polarons are created by the double-exchange mech- 
anism in order to gain kinetic energy, and if the AF 
couplings between the neighboring sites are J,Jd = 0, 
the ground state is FM. We also find evidence for a FM 
ground state for very small values of J, Jd < 0.05i for all 
models. The absence of interactions among the polarons 
in the symmetric model lead to a finite spin gap in that 
case. 

However, even for the asymmetric model the holes are 
only a weak perturbation of the underlying spin liquid. 
A hierarchy of energy scales is thus given in the spin sec- 
tor of the asymmetric model by the lowest-lying, gapless 
AF interactions among the polarons and second, larger 
energy scale of the order of the spin gap of the undoped 
system by the underlying spin liquid which remains in- 
tact. As a consequence of the spin liquid background, 
the spin-spin correlations decay exponentially fast for all 
models for the system sizes in our study (up to 2 x 256 
sites). However, we expect a crossover to power- law de- 
cay for the asymmetric at large distances. Also the string 
correlation function, which quickly approaches a nonzero 
value of the long-ranged order in the undoped system, re- 
mains finite for all models. Again, we expect a crossover 
to power-law decay for the asymmetric model at large 
distances. 

Since hole pairing on the rungs is only possible for the 
symmetric model without strong Coulomb repulsion, we 
find completely different dominant correlations functions 
in thermodynamic limit of infinite system size. The sym- 
metric model shows dominant pairing correlations and 
K p > 1 without Coulomb repulsion, in contrast to dom- 
inant CDW correlations with K p < 1 for the symmetric 
model with Coulomb repulsion and for the asymmetric 
model. In the FM region for small values of the AF in- 
teractions J, Jd < 0.05t, the triplet pairing correlations 
are strongly enhanced. 
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